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ABSTRACT 

In chis paper, Che probability density 
functions (pdfs) of the co- and cross-polar lzed 
phase differences are derived for backscatter from 
vegetation using the coherent and incoherent 
scattering theories. Unlike previous derivations, 
no assumptions or observations other than the 
applicability of the Central Limit Theorem (CLT) , 
the low fractional volume of the medium, the 
reciprocity of the scatterers, and the azimuthal 
symmetry of the scatterer’s orientation 
statistics, are employed. Everything else follows 
logically via the mathematics. The difference 
between the coherent theory and the incoherent 
theory is referred to as the backscatter 
enhancement effect. The influence of this 
enhancement effect on the phase difference pdfs is 
examined and found to be important under combined 
conditions of scatterer anisotropy and appropriate 
reflection coefficient values. 

Keywords: Polar ime trie , Backscatter, Vegetation, 
Phase Statistics. 

INTRODUCTION 

The phase difference statistics ■ of 
polarimetric radar returns from vegetation have 
been a subject of growing interest in the remote 
sensing community. This is due to the possibility 
that a strong dependence between the biophysical 
parameters and the measured phase difference 
statistics of the backscattered electric field 
components exists. Eom and Boemer [1] first 
derived the single-look co-polarized phase 
difference pdf as a function of one parameter: the 
amplitude of the correlation coefficient of the 
two scattered field components. This was 
subsequently generalized by Touzi and Lopes [2] to 
include a second parameter: the phase of the 
correlation coefficient. A multi-look 
co-polarized phase difference pdf was later 
derived by Lopes et al. [3]. Both co- and 
cross-polarized phase difference pdfs have been 
derived by Sarabandi [4] in terms of the Mueller 
matrix elements. Again, the co-polarized pdf is 
shown to be completely specified by two 
parameters: the degree of correlation and the 
polarized phase difference, as named in [4]. 

In this study, the co- and cross-polarized 
phase difference pdfs for backscatterlng from 
vegetation are derived using the coherent and 
incoherent theories, and the two are compared. 
The vegetation is modeled as a layer containing a 
random distribution of uncorrelated discrete 
scatterers over a flat surface. By considering a 


large number of scatterers , the CLT is employed 
and a multivariate Gaussian distribution for the 
real and imaginary parts of the backscattered 
fields results. The analysis confirms certain 
assumptions in [4] that appear to be made from 
observations; namely, the independence of the 
co-polarized and cross-polarized scattered fields. 
By assuming reciprocity of the scatterers and 
azimuthal symmetry of the scatter’s orientation 
statistics, we show the latter to be true for 
dipole scatterers and conjecture its correctness 
for larger scatterers. In addition, the two 
parameters which completely specify the 
co-polarized phase difference pdf are related to 
analytical expressions obtained via each theory. 

Since the coherent (or Distorted Born 
Approximation (DBA) ) theory is based on the 
addition of fields and the incoherent (or 
first-order Vector Transport (VT-1) ) theory is 
based on the addition of powers, certain 
interference terms appear in the DBA theory which 
do not appear in the VT-1 theory. These terms are 
due to the coherent interaction of 
counter-propagating fields following a 
scatterer-ground path, and give rise to the 
backscatter enhancement effect. Consequently, the 
results considered will look at how the 
enhancement effect influences the phase difference 
pdfs for model parameters which simulate 
vegetation canopies. In addition, since the 
co-polarized phase difference pdf is completely 
specified by two parameters, these parameters will 
also be examined as the angle of incidence, for 
example, is varied. 

PROBLEM FORMULATION 

The co- and cross-polarized phase difference 
pdfs for the elec tromagnecic backscatter from 
vegetation are derived in this section. The 
vegetation is modeled as a layer of thickness d 
consisting of a sparse distribution of identical 
electrically-thin lossy scatterers. The 
scatterers are assumed to be uncorrelated and 
uniformly distributed azimuthally, but can be 
assigned arbitrary elevation statistics. The 
underlying ground is represented by a lossy 
dielectric half-space having a flat surface. 

The resulting backscattered field with 
polarization p, due to an incident field of 
polarization q, can be expressed In the radiation 
zone as 

S' (s.l) - exp(ik|*|) E ( 4 ) , p,q € {h.v} (1) 
pq p<i 
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and h and v represent horizontal and vertical 
polarization, respectively. The variable, 

e (s ) , in (2) represents the scattered field 

pq r\ 

pattern for a single scatterer located at 5.^ in 

the layer containing N scatterers. Both E^(i) 

and e (i ) are complex quantities, 
pq « 

As in [ l ] - r 4 ] , che Central Limit Theorem 
( CLT) is applied. The CLT states that the pdf of 
a sum of independent random variables (rvs) 
approaches a Gaussian as N->». Hence, if assume 


that the X * s , as we 


11 as the Y 's, are 


independent rvs, then by the CLT, X^ and Y^ are 

each Gaussian rvs . Furthermore, since p and q can 
each assume two different states, a total of eight 
Gaussian rvs result. 

Assigning X, - X^, X^ - Y^, X^ - X^, X^ 

Y X-X X-Y.X-X and X - Y , 
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the icint Caussian pdf can be written as 
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where A is the covariance matrix and is given by 

A - (<X X >) , i, j e (1 3) ( 4 > 
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and | A | denotes its determinant. Also, x T denotes 
the transpose of the column vector x - ^ — 

Using (2), it is observed that 
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and a (d), a (dri), and a (drc) represent, 
respectively, che direct. direct-reflected 
incoherent, and direct-reflected coherent 
contributions to the total scattering coefficient. 
The three contributions are shown in Figure 1. It 
should be mentioned that (6) becomes of order * 
(i.e. o(«)) for very thin slabs and consequently a 
more detailed analysis, retaining terms of order 
e, is needed for this special case. 



Figure 1. The three different cor tribut ions to the 

o 

total scattering coefficient a 

To zeroth order in t, then, (5) with (6) and 
(7) becomes 
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Explicit analytical expressions have Q been 
derived for the different contributions to a for 
all polarization combinations, and can be found 
scattered in several papers (see Chauhan et. 
al.(6], for instance). Below, for the purpose of 
demonstrating the independence of Che co- and 
cross-polarized scattered fields, we give che 
analytical expression for the direct contribution 
to a (p>*q) ■ i ■ e • 
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However, by che DBA theory (which is valid 
for small fractional volume, e « 1), w « have 

the analytical results 

<E E * > - (A/4w) a° , + o(«) (6) 

pq p'q' p^p * 


<E E > - o(<) 
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a ° - (d) + o° (dri) + a° (drc) (8) 
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Here, a° is che total polarimetric scattering 
coefficient and is defined in the regular manner 
as (see, for example, Borgeaud, et. al.[5|) 
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Here f (-£’,£*) represents the scattering 

pq 

amplitude of a single scatterer with incident 
radiation of polarization q in a direction x - anc * 
scattered radiation of polarization p in a 
direction -£* (i.e. the backscatter direction); 

the overbar denotes the average with respect to 
orientation; the variable p represents the density 
of scatterers; and d represents the vegetation 
layer thickness. 

It can be shown, at least for dipole 
scatterers (but probably also for larger 

scatterers), that a (d) - 0 (p + q) • In a 

pppq 

similar manner, a (dri) - o (drc) - 0. In 
pppq pppq 

addition, by reciprocity, o - a - o - 



0 (p H q) also holds. Thus, the potarimecric 

matrix P which results from the DBA theory 

reduces, for backscatter, to 
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Making use of (10) and (13), the covariance 
matrix of (4) becomes 
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resulting co-polarized pdf, after integrating out 
p and p , is a function of only the 

w hh 

difference <£ - <p - <0 and can be written as 
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The parameters a and Q completely specify the 
pdf (18) and depend on the results obtained via 
either theory. The parameter a (0 < a < 1) 
denotes the degree of correlation between the 
co-polarized return signals and is related to the 
width of the pdf. The parameter & (-* 5 3 < *). 
on the other hand, has been called the polarized 
phase difference [4] and is related to the 
positioning of the center of the pdf (18). 


RESULTS 


and the joint Gaussian pdf can now be expressed as 
the product 
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where f^ and f^ are the joint Gaussian pdfs 

eo cross 

for the co-polar ized and cross-polarized rvs , 
respectively. It is evident front (15) that the 
two sets of rvs are independent. 

The cross-polarized covariance matrix given 
by (14b) is reducible to a 2x2 diagonal matrix 
with < 7 ° as the non-zero elements. This implies 
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where X and Y are independent rvs with equal 

vh »h 

variance . 

Making the polar coordinate transformations 


X - fi cos^ and Y - p sin^ (17a) 

vh vh vh vh vh vh 


X - p cos^ and Y - p stn^ , p € (h,v) 

PP PP 9V PP PP PP 

(17b) 

it is found that p^ is Rayleigh distributed and 
4 is uniformly distributed, the two being 

vh 

Independent. The conditional phase difference 
- $ ) pdf Is thus always uniform. The 
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The co-polarlzed phase difference pdf given 
by (18) is examined in this section for both the 
coherent (DBA) theory and the incoherent (VT-1) 
theory. Recall that the difference between both 
theories stems from the coherent direct-reflected 
contribution, < 7 °(drc) , caused by counter- 
propagating waves following a scatterer-ground 
path. It is this contribution which gives rise to 
the backscatter enhancement effect. 

The co-polarized phase difference pdf, via 
each theory, is shown in Figure 2 for a 
vegetation-like canopy consisting of a 
distribution of elliptical discs. The modeL 
parameters are given in Table 1 . Note from Figure 
2, that there Is a marked difference between the a 
parameters obtained using the coherent theory 
(i.e. a ) and the incoherent theory (i.e. a). 

e 1 

As mentioned earlier, a (given by (20a)) is 
one of two parameters needed to fully specify the 
pdf, and is related to the width of the pdf. It 
Is interesting to note, via (20a), that a ^ 

when either the direct contributions of all the 
< 7 °s are much greater than the respective 
direct-reflected contributions (i.e. |tf°(d) | » 
(<7°(dr)| ), or vice versa (i.e. |<7°(dr)| » 
|<y d (d)| ). Thus, in order for a + a^, a 

significant direct-reflected contribution must be 
necessary for at least one of the as in the 
denominator * of (20a). 

For the particular situation presented in 
Figure 2 with the parameters listed in Table 1, 
the elliptical disc scatterers are very long and 
thin, and are distributed with a near-hor izontaL 
elevation angle pdf. In addition, the dielectric 
constant of the ground is high. This means tha- 



che scattering which takes place due to the discs 
has a strong preference for horizontal 
polarization, and this, combined with a relatively 
high horizontally polarized reflection 
coefficient, results in a high direct-ref Lec ted 
contribution for cr^ (relative to the direct 

contribution c^( d^ and the total contribution 
) , and consequently a substantial difference 


between and a . Figure 3 shows the values of 

a and a as the angle of the incident wave 

impinging on the canopy varies. Note^ that the 
difference is about 0.2 for 55. This is 

twice the difference for 0-45° (the situation of 


Figure 2) . 

In conclusion, it can be said that the pdf 
obtained from the coherent theory may be very 
different from the pdf obtained from the 
incoherent theory; or, more simply, that the 
backscatter enhancement effect may have an 
important affect on the resulting phase 
statistics . 


_ (08*) THEORY : 3 , - 3*9. J, - -22.2* ’ 

^COHERENT (vr-1) THEORY a, * = -22.8* 

1 



Figure 2. The co-polarized phase difference pdf 
for backscatter ing from vegetation 
characterized by the parameters in Table 
1: DBA and VT-1 theories. 



Figure 3. The variation of the a parameter of the 
co-polarized phase difference pdf vs. 
angle of incidence: DBA and VT-1 

theories . 


Vegetation Model Parametersi! 

- ; 

Value 

General parameters: 


Frequency, f 

1 . 5 GHz 

1 Angle of incidence, ^ 

45° 

: 1 

Vegetation layer j 

thickness, d 


Dielectric constant of 

ground, t j 

t 

! 

1 

23.0+114.1 

* 

| i 

i 

Scatterer parameters: 


Elliptical discs with 


semi-major axis, a 

10 cm 

i 

semi-minor axis, b | 

1 cm 

thickness, t 

0.2 mm 

dielectric constant, « 

c 1 

17.0+15.55 

* 1 

density, p j 

4000/m 1 

elevation pdf, p(0) 

’ [ 4 sin 80 | 

j { 0 , otherwise i 

j | 

Table 1. Model parameters 
Figures 1 and 2. 

used to generate 
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